home *** CD-ROM | disk | FTP | other *** search
/ MacFormat 1994 August / August CD.bin / Shareware / Education / numericalmethods Folder / chap_7 / rctrap.m < prev    next >
Encoding:
Text File  |  1994-06-05  |  526 b   |  23 lines  |  [MATF/MATL]
  1. function T = rctrap(f,a,b,n)
  2. % T = rctrap(f,a,b,n)
  3. % Quadrature using the recursive trapezoidal rule.
  4. % f is the name of the function, input.
  5. % a is the left  endpoint, input.
  6. % b is the right endpoint, input.
  7. % n is the number of times for recursion, input.
  8. % T is the recursive trapezoidal rule list, output.
  9. m = 1;
  10. h = b - a;
  11. T = zeros(1,n+1);
  12. T(1) = h*(feval(f,a) + feval(f,b))/2;
  13. for j = 1:n,
  14.   m = 2*m;
  15.   h = h/2;
  16.   s = 0;
  17.   for k=1:m/2,
  18.     x = a + h*(2*k-1);
  19.     s = s + feval(f,x);
  20.   end
  21.   T(j+1) = T(j)/2 + h*s;
  22. end
  23.